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ABSTRACT 

Context. Dust attenuation in galaxies is poorly known, especially at high redshift. And yet the amount of dust attenuation is a key 
parameter to deduce accurate star formation rates from ultraviolet (UV) rest-frame measurements. The wavelength dependence of 
the dust attenuation is also of fundamental importance to interpret the observed spectral energy distributions (SED) and to derive 
photometric redshifts or physical properties of galaxies. 

Aims. We want to study dust attenuation at UV wavelengths at high redshift, where the UV is redshifted to the observed visible light 
wavelength range. In particular, we search for a UV bump and related implications for dust attenuation determinations. 
Methods. We use photometric data in the Chandra Deep Field South (CDFS), obtained in intermediate and broad band filters by the 
MUSYC project, to sample the UV rest-frame of 751 galaxies with 0.95 < z < 2.2. When available, infrared (IR) Herschel/PACS* 
data from the GOODS-Herschel project, coupled with Spitzer/MlPS measurements, are used to estimate the dust emission and to 
constrain dust attenuation. The SED of each source is fit using the CIGALE code. The amount of dust attenuation and the character- 
istics of the dust attenuation curve are obtained as outputs of the SED fitting process, together with other physical parameters linked 
to the star formation history. 

Results. The global amount of dust attenuation at UV wavelengths is found to increase with stellar mass and to decrease as UV 
luminosity increases. A UV bump at 2175 A is securely detected in 20% of the galaxies, and the mean amplitude of the bump for 
the sample is similar to that observed in the extinction curve of the LMC supershell region. This amplitude is found to be lower in 
galaxies with very high specific star formation rates, and 90% of the galaxies exhibiting a secure bump are at z < 1 .5. The attenuation 
curve is confirmed to be steeper than that of local starburst galaxies for 20% of the galaxies. The large dispersion found for these 
two parameters describing the attenuation law is likely to reflect a wide diversity of attenuation laws among galaxies. The relations 
between dust attenuation, IR-to-UV flux ratio, and the slope of the UV continuum are derived for the mean attenuation curve found 
for our sample. Deviations from the average trends are found to correlate with the age of the young stellar population and the shape 
of the attenuation curve. 
Conclusions. 

Key words, galaxies: high-redshift - galaxies: ISM-galaxies: starburst-ultraviolet: galaxies-dust: extinction 



1. Introduction 

Although dust is a minor component in galaxies by mass, its ef- 
fect on the observation of their stellar populations is striking. 
Dust captures a large fraction of the stellar emission, especially 
at short wavelengths. This process makes the direct observation 
of stellar populations from the UV to the near-IR, where they 
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* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



emit their light, insufficient to recover all the emitted photons. 
Reliable dust corrections are mandatory for measuring the star 
formation rate in the universe and its evolution with redshift 
from UV-optical surveys. When dust emission is measured, ac- 
curate star formation rates can be derived by combining IR and 
UV data, but these data are often not available, in particular for 
deep optical surveys (e.g., |Ilbert et al.| |2010). As a consequence, 
it is particularly important to study the dependence of dust atten- 
uation on parameters such as the observed luminosity, the stellar 
mass, or the slope of the UV continuum, since it could be used 
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to correct large samples for the effect of dust attenuation, at least 
in a statistical way. 

Any modelling of stellar populations in galaxies must also 
include attenuation from interstellar dust. Solving the radiation 
transfer in model galaxies is the best way to build physical and 
self-consistent SEDs. These models calculate the effective ob- 
scuration and produce attenuation curves as outputs, which are 
generally very different from the extinction curves affecting the 
flux of each star (e.g.JWitt & Gordon! [2000] IPierini et al.112004 



Tuffs et al. , 2004 Panuzzo et al. |2007| l. However, these sophis 
ticated models rely on numerous free parameters and physical 
assumptions that are difficult to constrain from the integrated 
emission from entire galaxies and for very large numbers of ob- 
jects. Simpler models have been specifically developed to anal- 
yse large samples of galaxies, introducing attenuation curves, 
recipes, and/or templates. The number of free parameters is rel- 
atively small. These codes are often developed to measure photo- 
metric redshifts and physical parameters such as the star forma- 
tion rate (SFR) and the stellar mass (M stel ). With the availability 
of mid and far-IR data for large samples of galaxies, new codes 
are emerging that combine stellar and dust emission on the basis 
of the balance between the stellar luminosity absorbed by dust 



and the corresponding luminosity re-emitted in the IR ( |da Cunha 
[eTaLl [20081 [Noll et al.||2009b) . 

Attenuation laws are introduced in all these codes, except 
for those which include a full radiation transfer treatment. The 



most popular attenuation curve is that of Calzetti et al. (2000 ), 
built for local starburst galaxies. Some specific recipes such as 
a time dependence of dust attenuation are sometimes introduced 
dCharlot & FalU |2000l [Panuzzo et al.||2007> . 

The attenuation law for local starbursts, based on spectro- 
scopic data, does not exhibit a bump at 2175 A such as that 
observed in the extinction curves of the Milky Way (MW) or 
the Large Magellanic Cloud (LMC) ( |Fitzpatrick & Massa[|2007| 
Gordon et al.||200'3j). The presence of a UV bump in the atten- 



uation curve of galaxies remains an open issue. In the nearby 
universe, the UV wavelength range has been investigated thanks 
to GALEX observations, but the results remain controversial. 



Wijesinghe et al. (GOTTb analyse the consistency of SFR indica- 



tors based on GALEX measurements in the far-ultraviolet (FUV) 
and near-ultraviolet (NUV) bands and fluxes in the H a line, and 
conclude that they must consider an obscuration curve without 
any 2175 A feature. Instead, from a careful analysis of pairs of 
galaxy SEDs, [Wild et al.| ( [20TT) conclude that the UV slope of 
the attenuation curve is consistent with th e presence of a bump 



at 2175 A, a conclusion also reached by Conroy et al. (2010 1 



from an analysis of the GALEX-SDSS colours of galaxies. At 
higher redshifts, several authors introduce a bump with mod- 



erate amplitude to improve photometric redshifts (Ilbert et al. 



2009 



Kriek et al. 201 1 1. Direct evidence of bumps comes from 

< z < 2.5: Noll and col- 
200 galaxies and 



the analysis of the galaxy spectra at 1 
laborators analyse high quality spectra of 
find significant bumps in at least 30% of the sources (Noll & 
|Pierini[ [2005] |Noll et al.| [2007] [2009a}. In an earlier paperTwe 



analysed SEDs of 30 galaxies in a redshift range between 0.95 
and 2.2, observed through intermediate band filters and with IR 
detections from HerschelfPACS, and found evidence for a UV 



bump in the dust attenuation curve of all these galaxies (Buat 



et al.||20T lb[ hereafter paper I). Selecting high-redshift galaxies 

immon (e.g. Lyman 
Daddi et ai.||2TX)5]), 



on their observed optical colours is quite common (e.g. Lyman 

_ |2008l~ 



break or BzK galaxies (Reddy et al. 



but most of the time mid- and far-IR data are not available for in- 
dividual targets, making it difficult to obtain any direct measure 



of dust attenuation. Studies of dust emission often rely on stack 
ing analyses (|Rigopoulou et al" 201 1| Burgarella et al. 2011 
Redd y et al.| |2012[ ). The sensitivity of Herschel in the deepest 
fields allows us to combine stellar and dust emission for individ- 



ual galaxies at high redshift ( Wuyts et al. 201 1 1, and to perform 
SED fitting accounting for both components (paper I). 

In this work, we will select galaxies at optical wavelengths 
with a very good sampling of their UV rest-frame SED: essen- 
tially, the sample is UV selected with a redshift range between 
1 and 2. When available, IR emission from Herschel/PACS will 
be added to better constrain dust attenuation. SED fitting will be 
performed to estimate the main characteristics of dust attenua- 
tion. The sample and the SED fitting process are described in 
Sects. [2] and [5] respectively. The global amount of dust attenu- 
ation and its variation with stellar mass and UV luminosity are 
discussed in Sect. [4] Section [B] is devoted to the description of 
the attenuation law. In Sect. [6f we revisit the relation between 
dust attenuation and the slope of the UV continuum. Our con- 
clusions are presented in Sect. [7] All magnitudes are given in 
the AB system. We assume that Q m = 0.3, £2a = 0.7, and 
H = 70kms-' Mpc 1 . 

2. Galaxy sample and rest-frame UV continuum 

2.1. Multi-wavelength data 



As part of the GOODS-Herschel key programme (E lbaz et al.| 
|2011| ), the Herschel Space Observatory ( |Pilbratt et al.[ |2010| ) 
surveyed part of the Great Observatories Origins Deep Survey 
Southern field (GOODS-S): 10' x 10' centred on the CDFS were 
observed at 100 and 160 /im over 264 hours by the PACS instru- 
ment (Poglitsc h et al.||2010| ). These are the deepest images of 
the universe ever obtained with PACS. 

Cardamone et al.| ( |2010| the MUSYC project) compiled a 



uniform catalogue of optical and IR photometry for sources in 
GOO DS-S, incorporat i ng the GOODS Spitzer IRAC and MIPS 
data ( Dickinson et al.[ [2003 ) and optical W FI broad-band data 



obtained at the ESO MPG 2.2 m telescope ( Hildebrandt et al. 
2006). In addition, they used deep intermediate-band imaging 
from the Subaru telescope to provide photometry with finer 
wavelength sampling than is possible from standard broad-band 
data. This was done to enable more accurate photometric red- 
shifts, and it provides a valuable means of tracing the detailed 
shape of the UV rest-frame spectrum. The GOODS-S field de- 
fined for the MUSYC project is larger than that observed by 
Herschel. In order to be able to put constraints on the emission 
of the sources in the far-IR, we restrict the field to that observed 
by the GOODS-He rschel programme. 

In paper I, we selected 30 galaxies that are strongly detected 
from the UV to the far-IR, and whose rest-frame UV was sam- 
pled by several intermediate-band filters. In the present work, 
we select a much larger sample of galaxies by relaxing the con- 
straint on the wavelength coverage and keeping galaxies without 
IR data. We start with the MUSYC catalogue, selecting sources 
with a spectroscopic redshift between 0.95 and 2.2. In this red- 
shift range we have more than ten photometric bands available 
in the UV rest-frame and a good sampling around 2175 A. As 
in paper I, we consider all the optical broad bands (U, f/38, B, 
V, R, I, z) and intermediate-band filters whose 5 cr depth was 
fainter than 25 ABmag (I A - All, 445, 484, 505, 527, 550, 574, 
598, 624, 651, 679, and 738 ( |Cardamone eTaT1 ( [20T0l Table 2); 
paper I, Fig. 1). We apply the same selection criteria as in pa- 
per I to ensure that sources have data below rest-frame 1800 A: 
a 5<x detection in the U and [738 filters for sources at z < 1.6, 
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1.4 1.6 
redshift 

Fig. 1: Redshift distribution of the galaxy sample. The distri- 
bution for the whole sample is plotted as a solid line, that for 
the galaxies detected by MIPS with a dashed line, and the filled 
histogram corresponds to the redshifts of sources detected with 
PACS. 



and in the B band and one intermediate filter corresponding to 
A < 1800 A in the galaxy rest-frame for z > 1.6. We further re- 
strict the sample to galaxies without a X-ray detection according 
to |Cardamone et aTlP008| l to minimise the contribution of an 
AGN to either the IR or UV emission of the sources. Our sample 
contains 751 sources, all detected by IRAC at 3.6 //m. 726, 398, 
and 353 sources are detected at 4.5, 5.8, and 8 /mi, respectively. 

We cross-correlate these 751 sources with the GOODS- 
Herschel catalogue obtained from source extraction on the 
PACS images performed at the prior position of Spitzer 24 /mi 
sources, as described in Elb az et aL] (j201 \\ . IRAC coordinates 
are used for the Spitzer 24 /mi sources, and we adopt a tolerance 
radius between IRAC and optical coordinates equal to 1 arcsec. 
290 sources have a 3 <x detection at 24 /mi. We check that all 
these sources have a single MIPS counterpart within a radius of 
3 arcsec. 76 sources are detected at both 24 /mi and 100 /mi, and 
30 at 24 /mi, 100 /mi, and 160 /an. 

During our analysis, we will distinguish between galaxies 
that do and do not have detections of dust emission at 24 /mi 
(hereafter MIPS data) and 100 /mi (hereafter PACS data). Given 
the small number of sources with 160 /an detections, we will not 
isolate the galaxies also detected at this wavelength from those 
detected at 100 /mi. 

The redshift distribution of the sources is shown in Fig. [T] 
The bimodality comes from the original distribution of spec- 
troscopic redshifts collected by Cardamon e et al.| ( |2010] > and 
mainly reflects the difficulty to measure spectroscopic redshifts 
for galaxies at 1.4 < z < 1.8. 55% of the galaxies selected have 
a redshift z < 1.5, and that percentage increases to 66 and 76% 
for the galaxies detected by MIPS and PACS respectively. 

2.2. UV continuum 



The UV continuum of star-forming galaxies is commonly mod- 
elled by a power-law, /i(erg crrT 2 s A ) oc where B is 
called the slope of the UV continuum. The observational deter- 
mination of B is made either from spectra or from broad-band 
photometry, sometimes with large uncertainties depending on 
the broad-band filters considered (Buat et al. , 201 lb). The wave- 




1.4 1.6 1.8 2.0 

log (F(24)/F(R)) 

Fig. 2: Fraction of L 1R attributed to an obscured AGN by 
CIGALE as a function of the flux ratio at 24 //m and in the R 
band, the error bars represent the standard deviation calculated 
by CIGALE for this parameter. The cross indicates the only ob- 
ject with F '(24) IF(R) > 200 and an AGN fraction higher than 
0.15 



length range over which B is measured also varies substantially 
( |Calzetti| [2001 ). Because of the very good sampling of the UV 
rest-frame continuum of our sources, we can accurately mea- 
sure the slope of the UV continuum from 1300 to 3000 A by 
performing a linear regression. The original wavelength range 
adopted by |Calzetti et al.|(|1994|l for fitting the UV continuum 
slope is 1200 - 2600 A. |Meurer et al. (1999} adopted a narrower 
range, 1200 - 1800 A. However, given our available photomet- 
ric data, extending the range to 3000 A adds extra photometric 
points, which enable us to better constrain the UV continuum 
slope without significantly biasing the results that would be ob- 
tained from a shorter wavelength range. We start at 1300 A to 
avoid the flattening of the UV continuum at shorter wavelen gths, 



down to 1200 A (Leitherer et al. 



2002 



Buat et al. 



2002). B is 



calculated in two ways, with and without data corresponding to 
the wavelength of the UV bump, and both estimates are found to 
be very similar. This is due to the fact that we have a large num- 
ber of photometric data for each source with only a few measures 
inside the bump area. In the following, we will use the values of 
B obtained with the full data set from 1300 to 3000 A. The stan- 
dard error on B is typically found to be of the order of 0.2. The 
rest-frame luminosity at 1530 A of each galaxy is obtained by 
interpolation, using the result of the linear regression and the 
spectroscopic redshift of the source. In the following, we will 
define Lfuv as vL v at 1530 A. 



(Code 



3. SED fitting with CIGALE 

3.1. The CIGALE code 

SED fitting is performed with the CIGALE code 
Investigating GALaxy Emission^] developed by Noll et al. 
(2009b). This is a physically-motivated code that derives prop- 
erties of galaxies by fitting their UV-to-IR SEDs. CIGALE com- 
bines a UV-optical stellar SED with a dust component emitting 
in the IR and fully conserves the energy balance between the dust 
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Table 1 : Input parameters for SED fitting with CIGALE 



Parameter 


Symbol 


Range 


Dust component 


distribution of dust mass 


(i- 


1, 1.5, 2, 2.5 1 


amplitude of the bump 




0, 0.5, 1, 1.5,2.0, 2.5,3,3.5 


full width half maximum of the bump profile 


7 


300, 350 A 2 


steepness of the attenuation law 


S 


-0.8, -0.6, -0.4, -0.2, 0, 0.2 


A v for the young stellar population 


A^(ySP) 


0.25, 0.5, 0.75, 0.9, 1.05, 1.2, 1.35, 1.5, 1.65, 1.8 mag 


reduction of A v for the old population 


/art 


0.25,0.5, l 3 


Stellar component 


age of the old stellar population 


t\ 


1,2, 3,4, 5Gyr 4 


e-folding rate of the old stellar population 


T\ 


l,3Gyr 


age of the young stellar population 


t2 


0.01,0.03, 0.1,0.3Gyr 


stellar mass fraction due to the young stellar population 


/ySP 


0.01,0.02, 0.05,0.1,0.2, 0.5 



Notes. Values of input parameters defining the dust emission and attenuation recipes (first part of the table) and the stellar populations (second 
part of the table). 

1 A single value of a (a = 2) is used when PACS data are not available. 

2 A single value of y is used for each run of the code. The baseline value is 350 A. 

3 A single value of / all is assumed for each run of the code. The baseline value is 0.5. 

4 The age of the old stellar population is always lower than that of the universe at the redshift of the object. 




Fig. 3: Comparison of parameters estimated by SED fitting for the catalogue of artificial galaxies. The initial values of the parameters 
are on the x-axis and the values estimated after SED fitting on the y-axis. The galaxies without IR data are plotted as black dots, 
the galaxies with only a MIPS detection as green filled circles, and those also detected with PACS as large red filled circles. The 
regression lines for each subsample are also plotted as solid lines and the 1 : 1 relation as a dashed line. From the left top to the right 
bottom, the parameters considered are the SFR, M star , L\ R , t\, logfe), log(/ y sp), A FV y, #b> and 6. 
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absorbed stellar emission and its re-emission in the IR. We re- 
fer to |Giovannoli et al.| ( |201 l| l for details on the use of CIGALE 
to fit SEDs of actively star-forming distant galaxies observed in 
their rest-frame UV and IR. In the present work, the sampling 
of the UV range is very good, so we expect to constrain the dust 
attenuation curve. To model the attenuation by dust, the code 
uses as a baseline the attenuation law of jCalzetti et al.| (|2000), 
and offers the possibility of varying the steepnes s of this law and 
adding a bump centred at 2175 A. We refer to|Noll et al](|2009b i 



for a complete description of the dust attenuation prescription. 
In brief, the dust attenuation is described as 



A(A) 



4.05 



(k'(A) + D Ao , y , Eb (A)) 



(1) 



where A v = 5500 A, k'(A) comes from |Calzetti et al.| ( |2000 1 
(Eq. 4) and Di , r ,E b (/l), the Lorentzian-like Drude profile com- 



monly used to describe the UV bump (Fitzpatrick & Massa 
[19901 |Noll et~aL)|2009al >, is defined as 



2„2 



D 



E b A z y 



(A 2 - Al) 2 + A 2 y 2 ' 



(2) 



The factor (j^j produces different slopes without modifying 
the visual attenuation Ay. It implies changes in the value of the 
effective total obscuration Ry, originally equal to 4.05 for the 
Calzetti et al. law ( |Noll et al.||2009b"l|Buat et al.||201 lb) . 

The peak amplitude above the continuum, full width at half 
maximum, and central wavelength of the bump (E b , y, and Aq) 
are free parameters of the coda 2 ] We fix the central wavelength 
of the bump Aq at 2175 A. The full width at half maximum of 
the feature y is not constrained by the SED fitting process. We 
adopt a value of 350 A intermediate between the value found 
for the MW and LMC and that found by |Noll et al.] ( |2009a) for 
high-redshift systems. We also ran CIGALE with y — 300 A and 
checked that the exact value of this parameter has no influence 
on the determination of the other parameters. The input values 
for the other parameters are listed in Table[T] The range of values 
for each free parameter was chosen after many trials: we began 
with a very large range of values and then reduced the range by 
excluding values never chosen during the y 2 minimisation. We 
will re-discuss the choice of the dust attenuation parameters in 
Sects. |]and[5] 

We adopt the stellar populations synthesis models of 
( 2005) . We consider a stellar population experienc- 
ing an exponentially decreasing star formation, whose age and 
e-folding time are free parameters. A younger stellar popula- 
tion with a constant SFR is added to the primary population. We 
adopt a Kroupa (2001 ) initial mass function (IMF). The two stel- 
lar components are linked by their mass fraction / y s P and an at- 
tenuation correction is applied to each of them. The same attenu- 
ation law is assumed for both stellar populations and a reduction 
factor of the visual attenuation (expressed in magnitude), / att , is 
applied to the old stellar population to account for the distribu- 
tions of stars of different ages ( |Charlot & Fall 2000; Panuzzo et] 
al. 2007) . For our baseline model, / att is fixed to 0.5, but we also 
performed runs with other values to check the validity of this 
assumption, since this parameter is ill-constrained (cf. Sect.|4|. 
Dust luminosities (Lir between 8 and lOOO/vm) are com- 



puted by fitting |Dale & Helou ( |2002) templates and are linked 



E b a nd y are related to the origi nal parameters c 3 and y FM intro- 



to the attenuated stellar population models (the stellar luminos- 
ity absorbed by the dust is re-emitted in the IR). The validity of 
Dale & Helou| ( [2002) templates for measuring total IR luminosi- 
ties of sources detected by Herschel is confirmed by the studies 
of |Elbaz et al!] l |20T0l|20rT) . When PACS data are not available, 
as is the case for the majority of the sources, only one of the tem- 
plates of |Dale & Helou| ( [2002) is used, which corresponds to the 
exponent of the distribution of dust mass over heating intensity, 
a, equal to 2. This is the mean value of a obtained for the sub- 
sample of galaxies with PACS detection. This template is very 
close to that used by |Wuyts et al. ( 2008) . Using a — 2 to extrap- 
olate the predicted fluxes in the PACS bands gives values consis- 
tent with the PACS flux limits from |Elbaz eTaT1 ( |20T7] l for 88% 
of the sample at 100 fim (flux limit: 0.8 mJy) and the percent- 
age increases to 92% at 160 /mi (flux limit: 2.4 mJy). Adopting 
different values of a modifies these percentages. With a - 2.5 
corresponding to more quiescent and colder galaxies, only 6% 
of the predicted fluxes for galaxies not detected with PACS are 
found above the flux limits at 100 and 160 /urn. For a hotter dust 
emission with a = 1.75, the fraction of galaxies which would 
have been detected by PACS reaches 18% at lOOyurn and 10% 
at 160/im. Therefore we do not favour low values of a for the 
galaxies not detected in IR. 

We run CIGALE on the full sample of 751 galaxies. The 
input parameters as well as additional output parameters like the 
stellar mass, the instantaneous SFR, or the total IR luminosity 
are estimated from the probability distribution function (PDF). 
The expectation value and its standard deviation are derived for 
each parameter. 84% of the sample is fitted with a minimum 
value of the reduced ;f 2 lower than 3. The list of free parameters 
is given in Table [T] 

Although we discarded sources with an X-ray detection 
(Sect. |2~T| ), highly obscured AGN could be present in our galaxy 
sample. To check this issue, we use the F(24)/F(R) flux ra- 
tio: highly obscured Compton thick AGN are expected to ex- 
hibit F '(24) /F(R) flux ratios larger than few hundreds (Fiore 



eTaT] [2009] l. We also run CIGALE adding an AGN contribu- 



duced by Fitzpatrick & Massa jl990) : E h = C3/yj; M and y = ypM X A\ 



tion as in paper I (with two obscured AGN templates from the 
Siebenmorgen et aQ ( |2004) library). In Fig. [2] we present the 
fraction of Lj R attributed to an obscured AGN by CIGALE as a 
function of F(24)/F(R). The highest value of F(24)/F(R) is 340. 
If we focus on galaxies with F(24)/F(R) > 200, only one source 
is found with an AGN fraction larger than 0.15. We discard this 
source from the following analysis. For the remaining sources, 
the average AGN fraction found by CIGALE is 0.09±0.07. Since 
this low value is consistent with no AGN contribution, we will 
use the results of the SED fitting process without any AGN tem- 
plate in the following analysis. 

3.2. Analysis of a catalogue of artificial galaxies 

In order to check the robustness and accuracy of the parameter 
estimation, we generate a catalogue of artificial SEDs from the 
input catalogue. First, we run the code on the full dataset to select 
the best model for each galaxy that corresponds to a single SED. 
The catalogue of artificial sources is built by integrating the best 
SED of each galaxy in the observed bands. A random error, typ- 
ically of the order of 10%, is added to the obtained fluxes. At 
the end, we have a sample of artificial galaxies whose properties 
and SEDs are expected to be representative of the real sources. 
The code is run on this artificial catalogue to compare the exact 
values of the physical parameters corresponding to the artificial 
SEDs to the parameters estimated by the code. With this cata- 
logue, we want to check our ability to estimate individual output 
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Table 2: Linear correlation coefficient (R) between the exact 
value of parameters and the value estimated by the code and its 
error (standard deviation cr of the PDF) of the estimated values 
for the sample of artificial galaxies. 



Parameters 


No IR data 


MIPS data 


PACS data 




R 


cr 


R 


cr 


R 


a 


\og(SFR) [M G yr'] 


0.86 


0.25 


0.99 


0.11 


0.99 


0.09 


log(M star ) [M G ] 


0.92 


0.21 


0.96 


0.15 


0.96 


0.15 


log(L m ) [L G ] 


0.84 


0.42 


0.98 


0.12 


0.99 


0.08 


h [Gyr] 


0.67 


0.98 


0.77 


1.02 


0.77 


1.00 
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parameters of the code given the wavelength coverage. We do 
not attempt to measure the effects of different signal-to-noise ra- 
tios. 

The results are summarised in Fig. [3] and Table [2] The 
sources detected with MIPS, with PACS and MIPS, and with- 
out any IR data are considered separately. For each parameter 
of interest in the present work, the correlation between the true 
and estimated values is shown in Fig. [3] The correlation coeffi- 
cient between both quantities is reported in Table [2] along with 
the typical uncertainty. All parameters are better estimated when 
IR data are included. No clear difference is found when PACS 
data are used in addition to MIPS data to estimate Lir: this is 
due to the way we choose the IR fluxes of the artificial galaxies. 
Indeed, the determination of Lir for artificial sources with MIPS 
data does not account for a variation of the IR templates, since 
artificial galaxies are built with only one IR SED (a = 2), also 
used to fit them. Therefore, Lir is well constrained with only one 
value. The SFR and Apuv strongly depend on the determination 
of Lir and are also well determined with only one IR (MIPS) 
data point. Without IR data, Lir appears to be slightly overes- 
timated, and as a consequence so are the SFR and Apuv- We 
will re-discuss the determination of Lir in the next subsection. 
The input parameters related to the star formation history are 
not well constrained: previous studies have also confirmed the 



weak constraints placed on these parameters ( Noll et al. 2009b 
|Giovannoli et al. 201 1| Buat et al. 201 la I. The parameters de 



scribing the dust attenuation curve, E\, and 6, appear better con- 
strained when IR data are available. The dynamical range of the 
estimated values is always lower than that of the true ones: the 
parameters are estimated as the mean of the PDF, which gives 
low weight to extreme values. In Sect. [5] we will re-discuss this 
effect, which strongly affects the determination of Lb. 

3.3. SED fitting with and without IR data 

It is important to check the robustness of the estimation of the 
most important parameters of this study against the availability 
of IR data. In complement to the analysis performed on artificial 
galaxies, we can also use the original sample with IR (MIPS or 
PACS data) to fit it with and without IR data and compare the es- 
timated values of these parameters. We consider the sources de- 
tected by MIPS (290 galaxies) adding PACS data when available 
(for 76 galaxies). We run CIGALE on the SEDs of these galax- 
ies without any IR data. In this case, only one IR SED is used 
with a = 2 and L lR is constrained by the obscuration of the UV- 
optical light. In Fig. [4] we compare the estimates of Apuv, Lir, 
Lb, and 6 with and without considering IR data. The agreement 
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Fig. 4: Parameters estimated with available IR data (x-axis) are 
compared to those estimated without any IR data. Galaxies with 
only MIPS data are plotted with small green symbols, galax- 
ies with MIPS and PACS detection with large red symbols. The 
green dashed line is the result of the linear regression and the 
black solid line is the 1:1 relation. Typical error bars (1 cr) are 
plotted on the left corner of each plot. 
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Fig. 5: Histogram of the distribution of the SFR obtained with 
SED fitting. The distribution with a dashed line corresponds to 
galaxies without IR detection, the solid line to galaxies detected 
by MIPS, and the filled histogram is obtained for galaxies de- 
tected by PACS. 



Fig. 6: Histogram of the dust attenuation in the FUV band. The 
distribution of A FUV is plotted with a dashed line for galaxies 
without IR detection, a solid line traces galaxies detected by 
MIPS, and the filled histogram indicates galaxies detected by 
PACS. 



between the estimates of all parameters is very good with corre- 
lation coefficients equal to 0.95 for Apuv, Lir, and E\,, and 0.92 
for 6. The SED fitting without IR data is found to slightly over- 
estimate dust attenuation and IR luminosities for small values of 
these parameters, by 0.27 mag for Apuv = 1 mag and by 0.15 dex 
for log(L IR ) = 10.5 [L Q ]. It is qualitatively consistent with the 
analysis performed on artificial galaxies with an overestimate of 
0.4 mag for Apuv = 1 mag and 0.20 dex for log(Ln0 = 10.5 [L Q ] 
for artificial galaxies without IR data (Fig. [3]). The discrepancy 
between exact and estimated values of Apuv and Lir is larger 
for artificial galaxies without IR data, with a slope between the 
true and estimated values of A FUV and L IR significantly different 
from 1 (see Fig. B), than between these quantities estimated on 
real galaxies observed in IR and analysed with or without IR data 
(see Fig.Q. Most of the pseudo-galaxies without any IR data are 
intrinsically faint when compared to the galaxies detected in IR 
and plotted in Fig. [4] This difference together with the nature of 
the data (observed and artificial fluxes) can explain the relative 
discrepancies found in the two analyses performed here and in 
Sect. |3.2| The determination of E\, is found to be robust against 
the inclusion or neglection of IR data, whereas 6 is more uncer- 
tain without IR data and is overestimated for low values of the 
parameter (by 0.1 units for 6 = -0.5). The latter result was also 
obtained with the catalogue of artificial galaxies. 

3.4. Star formation history 

As shown in Sect. |3.2| the detailed star formation history is 
not well constrained. The old stellar population is always found 
older than 1 Gyr, from ^1.5 Gyr at z — 2 to =s 3 Gyr at z — 1 . The 
average mass fraction locked in the young stars (/ y sp) is found to 
be of the order of 10%, which corresponds to galaxies with ac- 
tive star formation. In this work, we focus on the analysis of the 
UV part of the spectrum, which is dominated by the emission of 
the young stars. The average age of the young stellar component 
is 40 Myr: a steady state in the production of the UV continuum 
is not reached (Leitherer & Heckman, 1995 ). The distribution of 
the instantaneous SFR is shown in Fig. [5] 



4. The amount of dust attenuation 

Several output parameters of the code quantify the amount of 
dust attenuation. In addition to the V band attenuation of the 
young stellar population, the total amount of attenuation in the 
V and FUV bands (for both stellar populations) are also calcu- 
lated. The distribution of Apuv is shown in Fig. [6] The values 
range from 0.5 to 5 mag, the largest attenuations being found for 
galaxies detected with PACS. The mean dust attenuation is equal 
to 1.6, 2.2, and 2.6 mag for the galaxies without IR data, with 
MIPS data only, and with MIPS and PACS data, respectively. 
Note that dust attenuation is likely to be slightly overestimated 
for galaxies without any IR detection (Sect. [3~2| i. 

The relative amount of attenuation between the young and 
the old stellar population, expressed in magnitudes, is measured 
by f m . As mentioned above, / att is fixed to 0.5 for our baseline 
scenario, but we also performed runs with / att as a free param- 
eter allowed to take 3 values: 0.25, 0.5, and L This parameter 
is known to be ill-constrained ( |Noll et ak] |2009b[ |Buat et al) 
201 la). The value 1 is never chosen for the best-fit model and the 
mean value from the probability distribution function is found to 
be 0.52 ± 0.06. As a further check, we also ran CIGALE by also 
fixing / att at 0.25 and 1. All results presented in this work are 
found unaffected by the exact value of this parameter. In par- 
ticular, the distributions of Apuv, Lir, M stm -, E\» and 6 are not 
significantly modified when f m is varied. 



In Fig. [7] Apuv is plotted versus Lfuv and M star . There is a 
decrease of the upper envelope of the dust attenuation for large 
Lfuv in agreement with other findings based on UV selected 



samples at similar redshifts ( 


Burgarella et al. 2007 Buat et al. 


2009t|Burgarella et al.||201 1| 


Heinis et al.||2012). For a given UV 



luminosity, dust attenuation increases with stellar mass. Average 
values of Apuv are calculated per bin of UV luminosity and are 
also reported in Fig. [JJ in spite of a large dispersion, the mean 
value of A FUV is found to decrease when Lfuv increases from 
2.4 ± 0.8 mag for log(L FUV ) = 9.7 [L Q ] to 1.5 + 0.4 mag for 
log(Lpuv) = 10.7 [L Q ], which corresponds to a slope -0.9 in 
Fig. [7] The trend still holds when only galaxies not detected in 
IR are considered with a less steep decrease (slope reduced to 
-0.6). The results are not modified if we account for the over- 
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estimation of Apuv for galaxies without an IR detection (cf. 
Sects. 3.1 and 3.2 1. So the general variation seen in Fig. |7]can- 
not be attributed only to the sources detected in IR, with MIPS or 
PACS. An inverse trend was reported at z ~ 2 for Lyman break 
galaxies by Reddy and collaborators (Reddy et al. 2008, 2010) 
(but see also Burgarella et al. ( 2011[ l). Interesting" 
consider galaxies with z > 1 .8 to mimic a z 



,ly, if we only 
2 selection, our 

subsample contains galaxies more luminous than 1O 1O L and 
the distribution of Apuv as a function of Lfuv is consistent with 
a flat one (although we do not exclude a slight decrease when 
Lfuv increases). |Reddy et al.| ( |2010[ l also have galaxies with 
Lfuv > 10 10 L Q , which can explain the apparent discrepancy be- 
tween their analysis and the present work. The increase of the 
mean (and dispersion) of dust atten uation for galaxies f ainter 
than 10 10 L o in UV is confirmed by Heinis et al. ( 2012} from 
a stacking analysis of Herschel data in the COSMOS field. 

At higher redshift, from studies of the slope of the UV con- 
tinuum, Bouw ens et al.| ( |2009| l and |Lee et al.| ( |20PT| i also found 
a decrease of the dust attenuation with the absolute UV magni- 
tude. However, Finkelstein et al. ( 201 l| l and Dunlop et al. ( 2012| > 
found no significant variation of the UV slope with the absolute 
UV magnitude. 

Dust attenuation is also found to be correlated with the stellar 
mass as already reported from the local universe up to high red- 



shifts (Ma rtin et al.[ 2007; Iglesias -Paramo et al .||20071|Buat et 
aLl|2009HPannella et al .| |2009t|Garn & Best 



2010] [Finkelstein 



et al.| |2011[ |Sawicki| |2012) . |Martin et al.| pOO/J and |Bua~ 



al. ( 2009[ > studied UV selected samples up to z ~ 1. They re 
ported variations of A F uv as a function of the stellar mass or 
/f-band rest-frame luminosity. The mean relations they found 
are shown in Fig. |7| after some transformations. The relation of 
Martin et al. ( 2007 1 for z = 1 give n in Fig. |7] was obtai ned for 
a Salpeter IMF and the models of |Bruzual & Chariot] ( |2003 1. 
We translate it to a Kroupa IMF by reducing the stellar masses 
by 0.24 dex ( |Arnouts et at)[2007]l. The contr ibution of TP-AGB, 
included in the models of Maraston"] ( |2005|), is assumed to fur- 
ther reduce the stellar mass by 0.2 dex ( Muzzin et al.| 2009). 
Buat et al. (2009 1 proposed a relation between log(L IR /Lpuv), 
the rest-frame Zf-band, and the FUV luminosities, and did not 
find any trend with redshift. The relation reported in Fig. [7] cor- 
responds to a mean FUV luminosity (log(L FU v)) = 10.2 [L Q ] 
representative of our current sample. The /f-band luminosity is 
transformed in stellar mass using the calibration of Arnout s"et al.| 
( 2007 ) modified by Ilbert et al. ( 20 1 1 and A FU v is deduced from 
log(LiR/LFUv) following Buat et al. ( 201 la I. Although the slope 
of these correlations are consistent with what we find here, the 
average dust attenuation at a given stellar mass in these previous 
works is lower than our present estimate. We also compare our 
results to the relation obtained by Pan nella et al.| ( [2009] ) at z = 2 
by comparing SFRs from radio and UV measurements. We ap- 
ply the corrections described above to translate their results to 
a Kroupa IMF, also adding the contribution of TP-AGB. Their 
relation is much steeper than the mean trend we find, leading to 
larger Apuv at high M stal , where IR data are available (i.e. when 
our estimates are the most reliable). Anyway, it is difficult from 
this comparison to draw a firm conclusion about the evolution of 
dust attenuation with redshift at a given stellar mass. Different 
methods were used to derive stellar masses and dust attenuation 
factors, and sample selections are also different. A UV selection 
like that used in this work is biased against high dust attenua- 
tion and stellar mass which might induce a shallower relation 
than for a mass- or SFR-selected sample. Moreover it is very 
difficult to estimate stellar masses, especially in blue, low-mass 
galaxies (Ilbe rt et al.| |2010| l. Stellar masses and dust attenua- 
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Fig. 7: Dust attenuation in FUV versus observed FUV luminosity 
and stellar mass. Galaxies with no IR data are plotted as small 
dots, those with MIPS data with medium dots, and those with 
MIPS and PACS detection with large dots. Upper panel: Apuv 
versus log(L F uvX colour coded with log(M star ). The average val- 
ues of A F uv per bin of log(L FU v) are overplotted with their stan- 
dard dispersion, calculated by accounting for the uncertainty on 
the determination of Apuv- Lower panel: Apuv versus log(M star ), 
colour coded with log(L FU v)- The average values of A FU v per bin 
of log(M stal ) are overplotted with their standard dispersion (blue 
solid line). Relations from the literature and discussed in the text 
are overplotted as short dashed line (M artin et al.| ( |2007] l), dashed 
line ([Buat et al.| ( |2009] >) at z - 1, and dot-dashed line ( |Pannella| 
|etal.| ( |2009| l)"atz^ 2. 

tion must be measured with the same method at all redshifts and 
within samples selected in the same way in order to be com- 
pared. Nevertheless, we can conclude that there is a net increase 
of dust attenuation at larger stellar masses and that the amount 
of dust attenuation also varies with the (uncorrected) UV lumi- 
nosity. 

5. The attenuation curve 

5.1. Estimation of E h and 6 

The dust attenuation curve is described by the parameters E\> and 
5(Eqs.[T]and[2]i. 

The distributions of these parameters for the best models, 
estimated as the mean of the PDF, are plotted in Fig. [8] The ex- 
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Fig. 8: Distribution of parameters of the attenuation curve. 
Upper panel: distribution of £ b values obtained for the best 
models (filled histogram) and from the analysis of the PDF (solid 
line). Lower panel: distribution of 5 values obtained for the best 
models (filled histogram) and from the analysis of the PDF (solid 
line). 
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Fig. 9: Distribution of E^. Thin line: E\, estimated as the mean 
of the PDF. Thick line: £ b values corrected for edge effects 
(see text for details). The regressions used for the corrections 
are: £ b (SEDfitting) = 0.26£ b (model) + 1.08, £ b (SEDfitting) = 
0.47£ b (model) + 0.91, and £ b (SEDfitting) = 0.56£ b (model) + 
0.74 for the samples without IR detection, with MIPS detections 
only, and with MIPS and PACS detections,respectvely. 



treme values of 6 are not commonly chosen for the best model 
obtained by x 1 minimisation. However, it is not the case for £ b 
for which the null value is well represented in the best fits. By 
definition, only positive values of £ b are possible and the out- 
put value, given by the code as the mean of the PDF, will be 
overestimated when the exact value is close to 0. The trunca- 
tion of the PDF, when its mean value is close to the extreme 
input values of the parameter distribution, induces an over- or 
underestimation of the parameter (see Fig.[3]and Sect. 3.2 1. The 
amplitude of the effect depends on the exact shape of the PDF. 
|Noll et aL ( 2009b ) discussed the case of a Gaussian distribution 
with a given standard deviation: for a standard dispersion of 1 
unit, typical of our results for Z? b (see below and Fig. [TO]), all 



the values of £ b lower than 1.2 are affected and the correction 
ranges from 0.6 for £ b = to for £ b = 1.2. Another way to 
quantify the effect, which is less dependent on the shape of the 
PDF, is to use the correlations obtained for the artificial galax- 
ies between the exact and the measured values of £ b . In Fig. [3] 
we performed linear regressions, splitting the mock sample in 
objects without IR detections, those detected by MIPS, and ob- 
jects detected by PACS. The exact value of E\, (^(model)) was 
chosen as independent variable. The resulting regression formu- 
lae are applied to the observed sample to correct the values of 
E\, obtained from the PDF of the SED modelling. The results 
are plotted in Fig.[9]and the regression formulae are given in the 
figure caption. As expected, the distribution of the corrected val- 
ues is found to be broader. Galaxies with unphysically negative 
values for £ b represent 6% of the sample, whereas all negative 
values are consistent with a null value of Zs b at the one sigma 
level. The mean value for the entire sample remains unchanged 
(see below). 

To further check the parameter ranges and the reliability of 
the results, we performed an additional run with a very large 
range of values for E h and 6: -4 < E h < 4 and -0.9 < 5 < 0.9 
with a loose sampling to limit the computation time. Negative 
values of E\, are found for only 3% of the best models. The frac- 
tion reaches 16% for the mean value of the PDF, but except for 
one object, all these negative mean values are consistent with 
E\, = at a one sigma level. The distribution of 6 is very similar 
to that obtained with the parameters of Table [T] This check con- 
firms the results found with the baseline run and the reliability 
of the range of values found for £ b and 5. 

From an analysis of the mean and standard deviation of each 
PDF, we can measure the fraction of galaxies for which the bump 
is securely measured. Mean values and standard deviations for 
£ b and 6 are plotted in Fig. [10] We adopt the same criterion as 
in paper I to select sources with a secure detection of the bump 
as well as those with 6 confirmed to be negative, which corre- 
sponds to a probability of at least 95% for £ b > and 6 < 0. 
Assuming a Gaussian distribution with mean and standard de- 
viation given by the fit, it corresponds to E\, - 2cr £b > and 
6 + 2cr s < 0. Using this criterion, 20% of the whole sample is 
exhibiting a secure bump. The percentage reaches 40% (resp. 
47%) for MIPS (resp. PACS) detections. If we apply the correc- 
tions for the truncation of the PDF discussed above, the fraction 
of secure detections of the bump increases to 35% for the whole 
sample and to 54% (resp. 50%) for MIPS (resp. PACS) detec- 
tions. A similar percentage (20%) of the whole sample shows a 
secure negative 5. The percentage is 41% for galaxies detected 
with MIPS and 43% for those detected by MIPS and PACS. Note 
that only four SEDs are fitted with a 5 value found securely pos- 
itive (with 5 - 2o~5 > 0. It is worth noting that the error on the 
estimation of E\, and 6 is lower when the galaxies are detected 
with MIPS (or PACS). 

From the above discussion, it is clear that we must be cau- 
tious when reporting mean values for £ b or 5. The mean am- 
plitude of the bump for the whole sample is (£ b ) = 1.6 ± 0.4. 
Applying the corrections for the truncation at the edge of the pa- 
rameter range (see above) only increases the dispersion but not 
the mean value {{E\>) = 1.6 + 1.1). The mean amplitude is in- 
creased to (Eb) = 1.8 for galaxies detected in IR. These values 
are larger than those obtained in paper I. The difference comes 
from the range of initial values adopted for £ b and the estima- 
tion of the parameter as the mean of the PDF. In paper I, this 
range was < £ b < 2, based on the fact that values larger than 
2 were never taken for the best models. In the present work, we 
extend the range of values up to 3.5 (corresponding to the MW 



9 



V. Buat et al.: Dust attenuation in high redshift galaxies 



value) to minimise the impact of the truncation of the PDF, and 
because larger values than 2 are also obtained for some of the 
best models (Fig. 10 1. While the distribution of values based on 
X 1 minimisation are similar (same mean and median: E\, = 1 .5), 
the values estimated from the PDF are larger in the present case 
and were likely to be underestimated in paper I. We ran the code 
with the parameters defined in the present work for the sample 
of galaxies used in paper I and obtained (E\,) — 1.6 ± 0.3. 

The mean value (E b ) =1.6 corresponds to 45% of the bump 
amplitude of the MW extinction curve and is similar to that of 
the LMC supershell region (Fitzpatrick & Massa 2007] [Gordon| 
et al.| |2003] >. The large dispersion we find is likely due to the 



combination of uncertainties in the determination of the bump as 
well as due to intrinsic differences among galaxies and different 
environments within a galaxy. 

The mean value of S for the whole sample is (<5) = -0.27 + 
0.17. No significant difference is found when galaxies with and 
without IR data are considered separately. The value found in 
paper I was less negative: running CIGALE on the first sample 
with the present parameters gives (6) = -0.22 + 0.18 compared 
to (S) = -0.13 + 0.12 in paper I. The difference between both 
estimates indicates the uncertainty: modifying the range of Et, is 
likely to imply a difference in the determination of 6. The mean 
value found here ((6) = -0.27) corresponds to an increase in 
UV similar to that of the LMC extinction curve (Sect.|6~2]i. We 



do not find any trend between E\> and 5 (Fig. 1 1 



5.2. Differential slope and E b 

The imprint of the UV bump in attenuated spectra is expected to 
lead to a change of the UV slope across the UV bump. Following 
Noll & Pierini (2005), who defined a spectral index 734, we split 
the photometric data in two subsamples corresponding to A < 
2175 A and A > 2175 A. We restrict the redshift interval from 1 . 1 
to 1 .8 in order to have at least five data points in each subsample. 
We calculate the values of {3 corresponding to each wavelength 
range and subtract the two values with a linear fit: Af3 = j3(A < 
2175 A) - /3(A > 2175 A). The standard error on this differential 
slope, obtained by adding quadratically the errors on each slope, 
is of the order of 0.7, i.e. quite large. 

A clear trend is found in Fig. [12] all the galaxies with a se- 
cure bump detection exhibit a negative value of A(J3). In order 
to have reference values, A(J3) is also calculated on the SEDs of 
the best models using the same photometric bands as for the real 
data. Models without bump in the attenuation curve correspond 
to A(J3) - 0.5. This is consistent with Fig. 



12 if we account for 



the large error found for A(J3) and the overestimation of for 
low values of this parameter (Sect. 5.1 1. Given all these uncer- 



tainties, this comparison has to be taken with care. However, our 
criterion E\, - 2cr £b > to select a secure detection of the bump 
appears to be valid and even conservative, since it corresponds 
systematically to negative values of A(j3). 



5.3. Discussion 

We can summarise the results obtained on the attenuation curve 
for our galaxy sample: a large range of values is found for the 
amplitude of the bump and the slope of the curve. Since this dis- 
persion is not significantly reduced by neglecting galaxies with 
the most uncertain values (see Fig. [TO] , we can state that these 
variations probably reflect a real diversity of attenuation curves. 
These results imply that SED fitting codes must allow for varia- 
tions of the attenuation curve to fit the UV continuum of galax- 
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Fig. 10: Mean (x-axis) and standard deviation (y-axis) for the 
estimated values of E\, (upper panel) and 5 (lower panel). In the 
upper panel, galaxies with a secure (> 2 cr) detection of the bump 
are located to the right side of the solid line. In the lower panel, 
galaxies with an attenuation curve confirmed to be steeper than 
the Calzetti one are located to the left side of the solid line (see 
text for details). The symbols are the same as in Fig. [3] 



ies at all redshifts and retrieve reliable attenuations ( jllbert et aT] 
|2009] >. 

The presence of a bump is confirmed in at least 20% of the 
sources and up to 35% if we account statistically for the trun- 
cation of the PDF of this parameter. A steeper attenuation curve 
than that of |Calzetti et al.|(|2000[ ) is also robustly found in 20% of 
our galaxies. Noll et al. (2009a) confirm the presence of a bump 
in ~ 30% of their sources, which is consistent with our analysis. 
The amplitude of the bump is plotted versus the parameters de- 
duced from the SED fitting in Fig. 13 No correlation is found 



with SFR or Lir and very weak ones may exist with Apuv (corre- 
lation coefficient R = 0.25) or M star (R = 0.36). E\, is found to be 
lower for younger stellar populations, higher fraction of young 
stars, high specific SFRs (defined as SFR/M stal ) and higher red- 
shifts. To check the robustness of these trends against SED fitting 
uncertainties for galaxies not detected in IR, we applied correc- 
tions to log(SFR), log(M star ), log(Li R ), A FU v, and E\, based on 
the regression lines plotted in Fi g. [3] and f ound no significant 
modification (not shown here). |Noll et aX] ( |2009a] > also found 
that galaxies with evidence of a 2175 A feature host older stellar 
populations than galaxies lacking an evident bump. 
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Fig. 11: Slope of the attenuation curve versus amplitude of the 
bump: symbols are the same as in Fig. [3] Galaxies with a bump 
detected with a probability higher than 95% are overplotted with 
crosses and those with secure negative values of 5 at the same 
level of confidence are inscribed in a large circle (see text for 
details). 
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Fig. 12: Amplitude of the bump of the attenuation curve versus 
the observed variation of the slope of the UV continuum across 



the bump. The symbols are the same as in Fig. 1 1 



The decrease of the bump amplitude in galaxies forming 
stars actively could be related to a destruction of bump carriers in 
intense radiation fields or in supernovae shock fronts. Galaxies 
with high SFRs are recognised to be more compact than galaxies 
populating the main sequence defined by the SFR-M stal relation. 
The PAH molecules detected in the IR SEDs of the former galax- 
ies exhibit lower equivalent widths (El baz et aT)|201 1| and these 
molecules are often proposed to be at the origin of the UV bump 
(Xiang et al. 201 1] >. A deficiency in the production of bump car- 
riers can also be invoked to explain the decrease of in small 
galaxies, forming stars actively. Most bump carriers or the pre- 
products of them should be produced by carbon-rich AGB stars, 
which appear about 300 Myr after the onset of the star formation. 
Therefore, compact, low mass galaxies with a short but intense 
star formation history are not good places for the formation of a 
strong UV bump. 



Conversely, radiative transfer models for a clumpy medium, 
mixing of dust and stars, and dust-clearing effects for older stel- 
lar populations predict a decrease in the amplitude of the bump 
and the effective attenuation even if the extinction law exhibits 



a large bump (e.g. Witt & Gordon 



Panuzzo et al. 2007 Noll et al. 



Pierini et al. 2004 



2000 

2007). If galaxies with such 



attenuation properties significantly contribute to the sample, the 
fraction of galaxies with dust properties that allow for the for- 
mation of a strong UV bump should be distinctly higher than 
observed. Models assuming a clumpy distribution of stars and 
dust also predict a general flattening of the attenuation curve, 
although an age-selective attenuation acts in steepening the at- 
tenuation curve ( |Inoue et al.| |2006[ |Panuzzo et al.| |2007[ ). A 
patchy distribution of the dust is needed to obtain an attenua- 
tion curve as grey as that of Calzetti et al.| (2000]) st arting with 
a SMC-like extinction curve ( |Witt & Gordon) |2000l |Charlot &| 
|Fall| |2000]>. The mean attenuation curve found for our sample 



galaxies is steeper than the Calzetti et al. (2000) curve: if a sig- 
nificant flattening of the attenuation curve is at work, then the 
extinction curve must be very steep in our galaxies. The dust 
grain size distribution plays a major role in the general shape 
of the extinction curve and a steep increase in UV, which is at- 

is 



2005) 



tributed to a deficit of large grains (e.g. Cartledge et al. 
not unlikely in high-redshift galaxies forming stars very actively. 

On the other hand, our sample is characterised by UV bump 
strengths and attenuation curve slopes close to those typically 
observed for the extinction curve in the LMC supershell region. 
Noll et al. (2009a) show that the typical bump widths of lumi- 



nous star-forming galaxies at these redshifts also agree with the 
mean value of the LMC supershell region. Hence, the LMC su- 
pershell region may be a model region for high-redshift galaxies 
with significant UV bumps. In that case, the attenuation curves 
have to be similar to the extinction curves and this can also be 
achieved if the dust tends to show a screen-like distribution. This 
assumption is supported by the observation of strong interstellar 



absorption lines in galaxies with significant UV bump (Noll & 



Pierini 2005 1 Noll et al. 2009a i, which can be interpreted by 



a high covering fraction of young massive stars by dusty gas 
clouds. Screen-like dust topologies can be produced by large- 
scale outflows, which is a reasonable assumption for vigorously 
star-forming high-redshift galaxies. 

6. Dust attenuation and UV slope 

The slope ji of the UV continuum has been widely used to mea- 
sure dust attenuation in galaxies since the works of Calzetti, 
Meurer, and collaborators. They showed that the intrinsic value 
of this slope does not vary a lot as long as a high star formation 
activity is maintained in galaxies ( |Leitherer & Heckman| [1995; 



|Calzetti| 2001 1: any deviation from this value is assumed to be 
caused by dust attenuation. This was confirmed by the correla- 
tion between the observed slope f3 and the amount of dust at- 
tenuation measured by Lir/Lfuv found for local starbursts (e.g. 



is more complex for non-starbursting galaxies (e.g 




Burgarella 


etal. 2005| Boissier et al. 2007 Hao et al. 2011 




Boquien et 



al.||2012[ ~and references therein), this relation seems to hold for 



high-redshift systems forming stars intensively, such as Lyman 



break or BzK galaxies ( Daddi et al. 
|Magdis etal.l[20T0l ). 



2007| Reddy et al. 2008 



We have a reliable measure of f3 thanks to the intermediate 
bands included in the analysis. The presence of a bump could 
modify substantially the value of /? measured with broad-band 
data if they cover the bump area (Burgarell a et al.| |2005[ [Buat| 
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Fig. 14: Dust attenuation in FUV as a function of 
log(L IR /L FUV ) = IRX. The symbols are the same as in 
Fig. [3] The relation of |Overzier et al. ( 201 1] > for local starbursts 
is plotted with a solid line, while the relation proposed by Hao| 
et al.| ( 201 l| l for local star-forming galaxies is plotted with a 
black dotted line. 



|et al.| |201 lb| l. However, our measure of f3 is insensitive to the 
presence of the bump (cf. Sect. 2.2 1. We restrain the following 



study to galaxies detected in IR. Even if the determination of 
parameters related to dust attenuation is satisfactory without IR 
data (Sect.[3]l, the Li R /L FUV -/6 diagnostic is very sensitive to vari- 
ations of the dust attenuation and star formation histories. We 
prefer to restrict our study to the most secure cases. 



Fig. 15: Dust attenuation as a function of (3 (same symbols as in 
Fig. [3J. The solid line corresponds to the relation found in this 
work: Apuv = 1-70 x (J3 + 2.5). The dashed line corresponds 
to the relation of |Overzier et al.|(|2011[) (th eir inner relation). 

for 24 fjm selected 
(201 la) for nearby 



The relation found by [Murphy et al. (201 1 



galaxies at 0.66 < z < 2.66 and by Buat et a. 



galaxies selected at 140 /im are plotted as blue solid and dotted 
lines, respectively. 



6.1. A FUV versus log(L IR /L FUV ) 



We adopt the formalism of |Hao et al.| ( |201l) (first proposed by 
Meurer et al.| ( |1999| l) to relate A FUV to logC^R/Lpuy) (called 
IRX hereafter): 



A FUV = 2.51og(l+ flFUV xlO IRX ), 
where a F uv = 0.595 for local starbursts (Overzier et al 



(3) 



20 



20TTT > 
lOn 
,q.[3] 



and 0.46 for normal star- forming galaxies (Hao et al. 
Fig. [14J the results of our SED fitting are compared to 
The local starburst relation is found to be consistent with our 
data, whereas the relation valid for nearby star-forming galaxies 
underestimates the dust attenuation in our high-redshift galaxies. 



6.2. A FUV -j6 relation 

When the intrinsic UV spectrum is modelled by a power law of 
exponent /3, the relation between the attenuated and unattenuated 
fluxes, /(/I) and fo(A) (expressed per unit wavelength), and the 
corresponding slopes ft and /3q becomes 



\og(m)/f(A 2 )) logiMAO/MM)) 



log(/h/U 2 ) 



log(A l /A 2 ) 



(4) 



The wavelengths A\ and A 2 are taken in the UV range, where the 
modelling by a power law is valid. Then, dust attenuation at Ai 
can be written: 



2.51og(A 2 Mi) 
l-A(A 2 )/(A(A l ) 



x 08 - #>) = C(Ai , A 2 ) x OS - #,), (5) 



where A(A 2 )/A(Ai) depends on the dust attenuation curve 
(Eq. [TJ. In order to be consistent with our SED fitting results, 



we adopt A\ = 1530 A (corresponding to the FUV effective 
wavelength) and A 2 = 3000 A to calculate C(A.\,A.£). Using 
1200 and 2500A modifies the calculations by less than 1%. 
We found that 6 = -0.27 + 0.17 for the entire sample which 
corresponds to C{A\,A 2 ) - 1.70 (with a range of 1.46 - 1.97 
corresponding to the standard deviation for 5). For comparison, 
C(Ai,A 2 ) = 2.3 for the pure Calzetti law (6 = 0), 1.3 for the 
SMC and 1.8 for the LMC supershell using the extinction curves 
of Gordon et al. ( 2003) >. We can estimate f3o from Eq. [5] Using 
C(A U A 2 ) = 1.70, we find (J3 ) = -2.5 ± 0.3, while (J3 ) = -2.7 
(resp. -2.3) for C{A\, A 2 ) = 1.46 (resp. 1.97), consistent with the 
value expected for a young stellar population aged of ~ 40 Myr 
(see Sect. |3.4| >. The corresponding relation between /3 and A F uv 
(A FUV = 1.70x06+2.5)) is plotted in Fig.[l5]along with the data. 
We recall that /? is calculated on the entire UV range (Sect. |2j. 
The relation of Ove rzier et al.| ( [20TT] l for local starburst galaxies 
(we consider their inner relation, very close to the original one 
obtained by |Meurer et al. ( 1999| l) is also found to reproduce the 
general observed trend between ft and A F uv with a formal scatter 
for both relations reaching 0.5 dex. Their relation corresponds to 
C(A\,A 2 ) = 2.05 and fio = -2.19 with the formalism adopted 
in this section. Flatter and very dispersed relations are found in 
the literature (e.g. |Buat et al.[ |2011a| |Murphy et~aTj|2011| and 
references therein). The relations found by |Buat et aL (201 la 1 
and |Murphy et al.| ( |201 1) for IR-selected galaxies at low and 
high redshift are reported in Fig. 15 Although a steeper attenua- 
tion curve may flatten the relation (cf. Eq.|5]>, |B~uat et al. ( 201 la) > 
show that a flat distribution with a large dispersion in the Apuv- 
/3 diagram is consistent with an attenuation curve only slightly 
steeper than the Calzetti one and variations in the star formation 
history with a substantial contribution of old stars. 
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Fig. 13: Variation of the bump amplitude £b- Symbols are the same as in Fig. [3] with crosses overplotted on sources with a secure 
detection of the bump. From left to right and top to bottom: is plotted versus log(SFR), log(M star ), log(Li R ), the age of the old 
stellar population t\ (Gyr), the age of the young stellar population log(?2) (Gyr), the mass fraction of the young stellar population 
log(/ySp), ^fuv ( ma gX the specific SFR log(SFR/M star ), and the redshift z. 



6.3. IRX-ft relation 

We combine Eq.[3]for starbursts and Eq.[5]to deduce the IRX-/? 
relation: 



IRX = log(10' 



,0.4x1.703+2.5) 



l)-log(0.595) 



(6) 



The resulting curve is plotted with the data in Fig. 16 The rela- 
tion proposed by Overzier et al. (201 1 1 for local starbursts is also 
shown. The modification of the slope of the attenuation curve 
(steeper than the Calzetti one and similar to that of the LMC 

supershell) and of the intrinsic UV slope (from 2.2 to -2.5) 

modifies substantially the IRX-/? relation. The data are found to 
roughly follow both laws but a low value of /3o is necessary to fit 
the small values of slope and attenuation. 

The standard deviation from the Overzier et al. relation 
reaches 0.27 dex and is reduced to 0.23 dex when the relation 



derived in this work is used. It is similar to that found by Meurer 
et al.|(1999] l for local starburst galaxies. At z ~ 2, |Reddy et al. 
(2010 found that typical Lyman break galaxies roughly obey the 
local starburst IRX-/? relation with a formal scatter of 0.38 dex, 
slightly higher than our value. This relatively low dispersion is 



probably due to our UV selection and the high star formation 
activity of our galaxies, since the original law was first pro- 
posed for local UV-selected starburst galaxies. Subpopulations 
of dusty, IR-luminous galaxies at low and high redshift are 
known to deviate from the original starburst relation with a 



higher IRX for a given slope ft (Goldader et al. 2002 Chapman 


et al. 


2005 


|Howell et al.||2010||Papovich et al.| |2006|). |Penner 


etal. 


2012 


l select dust-obscured galaxies (DOGs) at z — 2 and 



use GOODS-Hrrschel data to study their 1R properties. They 
find that DOGs lie on the left side of the IRX-/? relation, with a 
systematic increase of IRX with ft. 

The influence of various parameters on the dispersion found 



in Fig. 16 is explored in Fig 17 There is a global increase of IRX 
with Lir for a given value of ft. Such a trend was already reported 
at low and high redshift (e.g. |Howell et al.| |2010| |Chapman et| 



aL] 2005 and references therein). We can also investigate the 
role of the variation of the attenuation curve (6) and of the star 
formation history (fto). As shown in Fig.[l7] the sources located 
above the mean law correspond to 6 values less negative than 
our mean law (S = -0.27): they are fitted with a greyer attenu- 
ation curve. Such a trend with 6 is easily explained: for a given 
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Fig. 16: The observed slope (3 is plotted versus log(LiR/Lpuv) = 
IRX. The symbols are the same as in Fig. [3] The solid line is 
the relation found in this work, while the dashed line is that of 
( [20TT) . 



Overzier et al. 



slope ft these galaxies experience a larger dust attenuation than 
predicted by our mean law and the fit gives a less efficient red- 
dening and thus a greyer attenuation. /3q is not an output parame- 
ter of CIGALE but its influence can be investigated by using the 
age of the young stellar populations as a proxy: galaxies with the 
lowest age (and therefore a small value of /?o) are located well 
above the mean law and conversely the oldest populations (cor- 
responding to large values of /?o) below it. As a consequence, 
the intrinsic value of /3o may also play a role in the locus of the 
galaxies in the IRX-/3 diagram (see also Boquien et al. (2012); 
|Heinis et al.|(|2012l)). 

Elbaz et al. ( |201 1\ relate the starburst activity to a higher 
specific SFR and a high compactness of the star formation. In 
Fig. [18] the distance to the mean IRX-/? law is plotted versus 
the specific SFR of our galaxies. A weak but positive correlation 
is found (correlation coefficient R = 0.20), which is stronger 
(R = 0.51) when only galaxies detected with PACS are consid- 
ered. This trend is consistent with a larger contribution of com- 
pact regions, which likely have greyer attenuation curves, with 



increasing specific SFR as proposed by Elbaz et al. (201 1 1. 



7. Conclusions 

We have studied the dust attenuation properties in a sample of 
751 galaxies with 0.95 < z < 2.2 observed with intermediate- 
band filters in order to sample correctly the UV rest-frame con- 
tinuum. IR data from MIPS (290 galaxies) and PACS (76 galax- 
ies) were used when available. Our results are: 

1 . SED fitting is performed with CIGALE. The dust attenuation 
law is described with two free parameters: its steepness (6) 
and the amplitude of a bump at 2175 A (Et,). The availabil- 
ity of IR data improves the determination of all the parame- 
ters related to dust attenuation. Global parameters such as the 
SFR, stellar mass, or the amount of dust attenuation are well 
determined, whereas the detailed star formation history is ill- 
constrained. Dust attenuations, total IR luminosities, and 6 
are found to be slightly overestimated in the regime of low 
attenuation and luminosity when IR data are missing. 

2. The dust attenuation in the FUV is found to have a large 
dispersion for low Lfuv (Lfuv < 10 10 L Q ). The mean value 





Fig. 17: The observed slope (3 is plotted versus log(Li R /L FUV ) = 
IRX. The dust luminosity Lir (top panel), 5 (median panel), and 
the age of the young stellar population (bottom panel) are 
colour coded. The solid line is the relation found in this work, 



while the dashed line is that of Overzier et al. ( 201 1 



of Apuv decreases from 2.4 to 1.5 mag when Lfuv increases 
from 10 9 7 to 10 10 ' 8 L Q . The dust attenuation is found to in- 
crease with stellar mass. 

The parameters describing the dust attenuation curve, E\, and 
6, are found to span a large range of values. This is likely 
to reflect intrinsic variations of the attenuation curve among 
galaxies and for different environments within a galaxy. The 
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-10.0 -9.5 -9.0 -8.5 -8.0 

log(SFR/M star ) (yr 1 ) 



Fig. 18: Distance of the data points in the IRX-/3 diagram from 
the mean law found in this work against the specific SFR. 
Symbols are the same as in Fig. [3] 



presence of a bump is confirmed in 20% of the sample, at 
Z < 1.5 for 90% of the confirmed detections. A dust attenu- 
ation law steeper than that of |Calzetti et al.| ( |20Q0] > (6 < 0) 
is also confirmed for 20% of the sample. These percentages 
increase to 40% for galaxies with IR detections. The mean 
values of £b and 6 derived for the whole sample are simi- 
lar to those found for the extinction curve of the LMC su- 
pershell. Finally, E\, is found to be anti-correlated with the 
specific SFR of galaxies. 
4. The relations between Apuv and /3 (the slope of the UV con- 
tinuum) and IRX (log(LiR/Lpuv)) and /3 are derived for our 
mean attenuation curve and for galaxies of our sample de- 
tected in IR. The intrinsic slope of the UV continuum, /3q, 
is found equal to -2.5 and the galaxies follow the IRX-/3 re- 
lation with a low dispersion. Galaxies found above the av- 
erage IRX-/? relation have a high IR luminosity and a dust 
attenuation curve slightly greyer than our mean law, which 
can be due to a higher compactness of the star formation 
regions. The age of the young stellar population (i.e. the in- 
trinsic shape of the UV continuum) is likely to also play a 
role: for a given (3, IRX increases when the age of the young 
stellar population decreases. 
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